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Abstract 



We discuss the measured partial width of the pseudoscalar charmo- 
nium state r]c into two photons. Predictions from potential models are 
CN I examined and compared with experimental values. Including radia- 

tive corrections, it is found that present measurements are compatible 
, both with a QCD type potential and with a static Coulomb potential. 

I The latter is then used to give an estimate on the rjf, decay into two 

• photons. Results for ijc are also compared with those from J/tlj data 

C^i through the NRQCD model. 

Oh! 

: 1 Introduction 

> . 

^ ■ In this paper we revisit the calculation of the two photon width of rjc, high- 

■ lighting newest experimental results and updating the potential model cal- 

culation. This allows for a reliable estimate of the two photon width of r]b, 
which is been searched for in 77 collisions [jl| . We shall see that the expected 
two photon width of r]b is within reach of the precision in the LEP data being 
analyzed. 

The charmonium spectrum has been the basic testing grounds for a vari- 
ety of models for the interquark potential, ever since the discovery of the J/i/j 
in 1974 [0. The experimental scenario describing the cc bound states is close 
to completion, with the observed higher excitation states ^Pq, ^Pi and spin 
2 states : decay widths into various leptonic and hadronic states have 
been measured and compared with potential models Most of this note 

is dedicated to examine the theoretical predictions for the electromagnetic 
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decay of the simplest and lowest lying of all the charmonium states, i.e. the 
pseudoscalar 77^. In Sect. 2 we shall compare the two photon decay width 
with the leptonic width of the J/ip., which has been measured with higher 
precision and found to be 15% higher than in previous measurements 0. 
This implies that a number of potential models whose parameters had been 
determined by the leptonic width of the J/ip may need some updating, and 
so do some predictions from these models. Potential model predictions for 
Vc 77 can be found in Sect. 3, together with a value for the two photon 
width of rji) extracted from the Coulombic potential. In Sect. 4 we show 
the predictions for r]c decay widths, using the procedure introduced in for 
the description of mesons made out of two non relativistic heavy quarks, by 
means of the Non Relativistic Quantum Chromodynamics-NRQCD. In Sect. 
5 we compare these different determinations with the experimental value of 
the rjc —>■ 77 decay width expanding some recent theoretical analyses on this 
subject (see for instance 0, and references therein). 



2 Experimental values and relation to J/ip 
electromagnetic width 

The first evidence of the rjc state has been found in the inclusive photon spec- 
tra of the ip' and J/ip decays [|T^. Subsequently, through 77 collisions, 
the decay width of rjc into two photons has been measured in different exper- 
iments. The most recently reported values for the radiative decay width are 
shown in fig. (|l|) |[T3|- p^] , together with the Particle Data Group average p3|1 , 
which reads 

rex.p(r?c^77) = 7.4±1.4keV. (1) 

In order to compare the experimental determinations with theoretical expec- 
tations, we start with the two photon decay width of a pseudoscalar quark- 
antiquark bound state |Q with first order QCD corrections which can 



be written as 

r(r/c ^ 77) = 



1 + -^ 



TT 

In eq. P^ is the Born decay width for a non-relativistic bound state which 
can be calculated from potential models. A first theoretical estimate for this 
decay width can be obtained by comparing eq. (0) with the expressions for 
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Figure 1: Experimental values of various measures of rjc 77 



the vector state J/ip [^], i.e. 

r(j/^ - 



e e 



B 



16 as 



(3) 



The expressions in eqs. (0) and @ can be used to estimate the radiative 
width of rjc from the measured values of the leptonic decay width of J/ip, if 
one assumes the same value for the wave function at the origin ip{0), for both 
the pseudoscalar and the vector state. This is true up to errors of 0{as/ml) 
(see for instance p7|,p^). 

Taking the ratio between eqs. and (0) and expanding in as, we obtain 



r(r/c ^ 77) 4 (1 - 3.38«,/7r) 



1 + 1.96— + 

IT 



(4) 



T{J/ip e+e^) 3(1- 5.34a,/7r) 
The correction can be computed from the two loop expression for as and 



the value |^ as{Mz) = 0.118 ± 0.003. Using the renormalization group 
equation to evaluate as{Q = Irric = 3.0 GeV) = 0.25 ± 0.01, and the latest 
measurement 

^expiJ/^ e+e") = 5.26 ± 0.37 keV (5) 

one obtains 



^{Vc 77) ± Ar(r/c ^ 77) = 8.18 ± 0.57 ± 0.04 keV 



(6) 
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where the first error comes from the uncertainty on the J/ip experimental 
width, the second error from . This estimate agrees within la with the 
value given in eq. (|1|). Here we assumed the scale to be Q = 2mc = 
3.0 GeV. This choice is by no way unique, and in fig. we show the 
dependence of the rjc photonic width, evaluated from eq. (H), upon different 
values of the scale chosen for a.. 
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Figure 2: The dependence of the rjc decay width to 77 (in keV) is shown 
with respect to the scale chosen for in the radiative corrections. The 
horizontal lines represent the central experimental value (full line) and the 
(dashed lines) indetermination. 



As one can see from fig. (|]) the experimental width value is not sufficient 
to uniquely determine the a scale choice of . We shall therefore include 
this fluctuation in the indetermination due to radiative corrections. 



3 Potential models predictions for r]c cind 775 

We present now the results one can obtain for the absolute width, through the 
extraction of the wave function at the origin from potential models. For the 
calculation of the wavefunction we have used four different models, namely 
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the Cornell type potential |^ V{r 



■ - + ^ with parameters a = 2.34 and 



k = 0.52, the Richardson potential |^ Vr 



r = --. 



127r 



with A^/ = 3 and A = 398 MeV, 



3 33-2Nf J (27r)^g2 log(l+q2/A2-) 

and the QCD inspired potential Vj of Igi- 



Ono 13T|,|2 



Vj{r) = VARir) + + ar, \/Ai?(r) = - 



4«(2)(r) 



(7) 



with two different parameter sets, corresponding to Aj^^r = 0.5 GeV and 
^Ms ~ 0.3 GeV^ respectively |]^. We also show the results from a Coulom- 
bic type potential with the QCD coupling frozen to a value of r which 
corresponds to the Bohr radius of the quarkonium system, obtained by solv- 
ing the equation tb = ?>/ {2mcQ.s{rB)) ||3^. We stress that the scale of 
occurring in the radiative corrections and the one of the Coulombic potential 
are different 
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Figure 3: The dependence of rjc decay width to 77 in keV for different 
potential models is shown function of rric . 



We show in fig. @ the predictions for the decay width from these po- 
tential models with the correction from eq. (0) at an scale Q = 2mc , 
observing that the calculated widths stay well within one standard deviation 
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of the width value given by eq. (|T]). For any given model, sources of error 
in this calculation arise from the choice of scale in the radiative correction 
factor and the choice of the parameters. Including the fluctuations of the 
results given by the different models, we can estimate a range of values for 
the potential model predictions for the radiative decay width T{ric — 77), 
namely 



r(r/c ^ 77) = 7.6 ± 1.5 keV . (8) 
ALEPH has recently started to search for r]b into two photons and it 



is interesting to see whether the potential models can predict for this decay a 
value within experimental reach. Predictions of course will be affected by the 
error due to the parametric dependence of the given potential model, an error 
which can be quite large since most of the parameters have been tuned with 
the charmonium system. On the other hand, the Coulombic potential gives 
results for the charmonium system in agreement with all the other models, 
and, at the same time, is relatively free of such parameter dependence. With 
only the scale of in the wave function to worry about, it can be used for 
a reliable estimate. For this purpose we shall make use of the expression in 
eq. (0) where this time Cq = 1/3 and mi, = 5.0 GeV. This gives the potential 
model prediction 

T{rib 77) = 0.50 ± 0.03 keV , (9) 

where the error is associated to different choices of rrib values and to the 
indetermination on occurring in the radiative correction. A check of this 
estimate can be given using the leptonic width of the T and the expansion 
given in eq. (H). To first order in one obtains: 



Tirik ^ 77) 1 (1 - 3.38a,/7r) 1 



, s - . , . 1 + 1.96— + (10) 

r(T^e+e-) 3 (1 - 5.34a,/7r) 3L vr ^'^J ^ ^ 

which differs from eq. (^ only by a charge factor; using the PDG average ||2^ 



re^p(T ^ e+e") = 1.32 ± 0.05 keV (11) 

and assuming the wavefunctions of the two bb bound states to be equal we 
have 

T{7]b 77) = 0.49 ± 0.04 keV (12) 

in agreement with the Coulombic model prediction eq. (^). For the radiative 
correction factor we have used as{Q = ^rrih = 10 GeV) = 0.18 ± 0.01. The 
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associated error in eq. (^) takes into account the indetermination on the 
experimental value eq. (|lTD and the one on • 

A more thorough discussion on the theoretical estimates of the rjf, de- 
cay into two photons and comparison to some recent calculations (see for 
instance [^) will follow in a future publication 



4 Octet component model 

We will present now another model which admits other components to the 
meson decay beyond the one from the colour singlet picture (Bodwin, Braaten 
and Lepage) 0. NRQCD has been used to separate the short distance scale 
of annihilation from the nonperturbative contributions of long distance scale. 
This model has been successfully used to explain the larger than expected 
J/ip production at the Tevatron and LEP. According to BBL, in the octet 
model for quarkonium, the decay widths of charmonium states are given by: 



T{J/ij LH) 



a, (1-3.7^ 

TT 



2{J/i;\O,i'S,)\J/ij)10al{n'~9) 



mi 

+ 3.2a (1-6.7— 

TT 



243 



X 



2(J/^|Pi(35i)|J/V^) 2 



a; (2.0a + 0.6) (13) 



V{JI^ e+e") 



87ra^(J/V>|Oi(=^^i)|J/V>) 
27m2 

327ra^(J/V;|Pi(^gi)|J/^) 
81m4 



16 a, 
3"^ 



i_ 



(14) 



^ LE) 
a 



X 



TT 



27rag(r/,|Oi(igo)|r/e) 

87ra'M\PiCSo)\ri,) 
27mf. 



/143 31 . 
1 + — - tttVT^ 



V 6 



24 



(15) 



r(r/c ^ 77) 



327ra2(r/,|Oi(i5'o)|r/c) 



81m^ 



I 4 / 3 TT 
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1287ra'{v,\P^{'So)\v,) 
243< ^'^^ 

There are four unknown long distance coefficients, which can be reduced to 
two by means of the vacuum saturation approximation: 

Gi = {JMO,('S,)\J/^P) = {vc\0,CSo)\Vc) (17) 

F, ^ {J/i,\P,{'S{)\J/i,) = {v,\P,CSo)\Vc) (18) 

correct up to O(f^), where v is the quark velocity inside the meson. We use 
the J/ip experimental decay widths as input in order to determine the long 
distance coefficients Gi and Fi . This result in turn is used to compute the 
Tjc decay widths. 

The BBL model gives the following decay widths of the rjc meson: 

r(r?c ^ 77) = 9.02 ±0.65 ±0.14 keV (19) 

and 

T{ric LH) = 14.38 ± 1.07 ± 1.43 MeV (20) 

where the ffist error comes from the uncertainty on the J/^p experimental 
width, the second error from ag ■ This results agree with experimental data 
within la, con&ming the applicability of the BBL model to the charm sys- 
tem. We leave to a future publication the application of the BBL to the 

6-system. 



5 Comparison between models 

For comparison we present in fig. (^) a set of predictions coming from different 
methods. We see that the theoretical results are in good agreement which 
each other. 

Starting with potential models, we see that the results are in excellent 
agreement with the experimental world average taken from PDG. The ad- 
vantage of this method is that we are giving a prediction from ffist principles, 
without using any experimental input. The second evaluation, given by BBL 
using the experimental values of the J/ ijj decay, is off by la from the central 
value. This is true also for the determination of the BBL model with non- 
perturbative long distance terms taken from from the lattice calculation ||38|| , 
affected from a large error. The advantage of the latter is that its prediction. 
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Figure 4: The rjc two photon width as calculated in this paper using (start- 
ing from below) Potential Models results, BBL model with input from J/ijj 
decay data, Lattice evaluation of Gi and Fi factors. Singlet picture with Gi 
obtained from J/ip ^ e^e~ and J/ip ^ LH processes respectively. The 
scale has been enhanced with respect to fig. 2 to evidentiate the difference 
between various theoretical results. 



like the one from potential models, does not make use of any experimental 
value. Next is the point given by the singlet picture from the electromagnetic 
decay of the J/ip, in agreement with the central value of the t]c. The last 
point is obtained also from the singlet picture with the J/ ip decay into light 
hadrons, and is in disagreement with the experimental measure. This is a 
long standing problem with some charmonium decay widths that hasn't been 
resolved yet (see for instance 13^ and references therein). 



6 Conclusions 

The T{r]c 77) decay width prediction of the potential models considered 
gives the value 7.6 ±1.5 keV which is consistent with the individual mea- 
surements and the world average p3 . The Coulombic model is in agreement 



with predictions from other models, and gives for the rih 77 decay width 
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the estimate 0.50 ± 0.03 keV. Predictions of the BBL model for the 7]^ 77 
decay width are consistent with the experimental measurements, for both the 
long distance terms Gi and Fi extracted from the J/ip experimental decay 
widths and the one evaluated from lattice calculations. 
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